Crossed Andreev reflection in a graphene bipolar transistor 
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We investigate the crossed Andreev reflections between two graphene leads connected by a narrow 
superconductor. When the leads are respectively of the n-and p- type, we find that electron elastic 
cotunneling and local Andreev reflection are both eliminated even in the absence of any valley-isospin 
or spin polarizations. We further predict oscillations of both diagonal and cross conductances as a 
function of the distance between the graphene-superconductor interfaces. 



Several decades after Einstein, Podolsky and Rosen 
raised their famous paradox the successfuU imple- 
mentation and study of polarization-entangled states of 
photons 0] has ruled out the possibility of simple local 
hidden- variables formulations of quantum physics [3] . In 
solid state physics, the controlled production and detec- 
tion of charge- or spin-entangled electronic states remains 
a major challenge, regarding the fundamental concepts 
of quantum physics, as well as quantum processing and 
communication issues. Owing to the structure of their 
ground state, conventional singlet superconductors were 
suggested as natural sources of spin-entangled 0, [E Si or 
even momentum-entangled electrons Q- Unfortunately, 
superconductors are also bad beam splitters since the 
electron-hole Andreev conversion is essentially a retrore- 
flection in usual metals or semiconductors [§| . Strikingly 
Beenakker uncovered that Andreev reflection (AR) may 
be specular in graphene 0, ■ Therefore it should be 
possible to observe paired electrons along diverging tra- 
jectories within a single graphene flake connected to a 
large superconducting electrode. Nevertheless angular 
filtering is a rather difficult task in quantum electronics 
in contrast to optics. Accordingly a lot of theoretical 
11, 12, [3] and experimental ^14] efforts have been de- 



voted to the crossed Andreev reflection (CAR) process 
by which a superconducting condensate (S) emits two 
quasiparticles in two normal metallic leads Ni and N2 
where they can be probed separately. The main draw- 
back of such N1SN2 junctions was identified as the ubiq- 
uitous presence elastic cotunneling (EC) and local AR. 
Indeed during the EC process an electron tunnels elasti- 
cally from Ni to N2 through the superconductor without 
any Cooper pair transfer, while in AR the paired elec- 
trons are injected in the same lead. In standard nonrel- 
ativistic conductors with low transparency tunnel con- 
tacts, the cross conductances originating from CAR and 
EC cancel exactly each other in the noninteracting limit 
ll( . and it is necessary to consider the noise properties 



to probe the CAR process |13j . 

In this Letter, we show that the unique relativis- 
tic band structure of graphene enables to observe a 
pure crossed Andreev reflection in a three-terminal n 
graphene/superconductor/p graphene (G1SG2) bipolar 



transistor, see Fig. [TJ Accordingly the injected Cooper 
pair is splitted in electrons which further propagate in 
opposite directions within Gi and G2 respectively. In- 
deed both EC and local AR may be totally suppressed 
owing to the presence of Dirac points in the spectrum 
of Gi and G2. In contrast to the nonrelativistic case, a 
CAR dominated transport should be observed directly in 
the conductance measurements performed on such bipo- 
lar graphene transistor (see Fig. I2I3I) w ithout resorting to 
noise [13| or interaction effects [15|, |l^| . Similar phenom- 
ena in usual conductors are prohibited by the fact that 
the corresponding Fermi energies are always much larger 
than the superconducting gap . By studying the inter- 
play of superconductivity [17[ with the very special dy- 
namics of massless relativistic quasiparticles at a bipolar 
pn junction [H, [H, 0, [2l|, [H, [2|| , we obtain the oscilla- 
tory behavior of both diagonal and cross conductances of 
the G1SG2 transistor as a function of the superconductor 
width. 

We consider a graphene sheet occupying the xy plane. 
A superconducting top electrode covers the region from 
X — to X — d, creating a proximity induced su- 
perconducting barrier (S) between the normal leads Gi 
{x < 0) and G2 {x > d). Moreover it was argued re- 
cently that metal coating might also induce supercon- 
ductivity in graphene [2J|. Due to valley and spin de- 
generacy, one may use a four-dimensional version of the 
Dirac-Bogoliubov-de Gennes equation [loj 
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functions {'^a+:'^b+) relative to one valley (-I-) and 
their time-reversed hole states attached 
to the other valley (— ). The indices A and B label 
the two sublattices of the honeycomb structure of 
carbon atoms. The kinetic Hamiltonian is given by 
VFP-(y — —iHvpicFxdx + ^ydy) where the Pauli matrices 
Ux and Uy act in the sublattice space as well as the iden- 
tity (Jo- The energy e is measured from the Fermi level 
of the superconductor and is the energy-independent 
Fermi velocity. The electrostatic potential U{y) in leads 
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FIG. 1: Top: Graphene-superconductor-graphene (G1SG2) 
transistor. We assume that a positive bias Vi is applied to 
Gi while S et G2 are grounded. Bottom: Incident electron 
at energy e = /i in n- doped graphene (Gi). The Andreev 
reflected hole (o) in Gi and the transmitted electron (•) in p- 
type G2 are "blocked" at the Dirac points since k' — 0. Thus 
the incoming electron may only be reflected as an electron 
(•) in Gi or transmitted as a hole (o) in G2 for any incidence 
angle a [25l |. For e ^ fi (not shown), the elastic cotunneling 
and the local Andreev reflection are still blocked provided a 
exceeds the critical angle ac(e) — arcsin(l(/x — e)| /{fj, + e)). 



Gi (i = 1, 2) and in the central region may be adjusted 
separately using state-of-the art local gates technology 
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22l |. It is assumed that U{r) = —fii and A(r) = 
in Gi, while U{r) = —^s and A(r) = Aqb"^ is finite for 
< X < d. This square-well model is fully justified by 
the unusually large Fermi wavelengths in graphene leads, 
and the fact that the Fermi wavelength beneath the 
superconductor should be far smaller, namely <C /ig. 

In order to clarify the physics of such bipolar G1SG2 
planar heterojunctions, we first give a simple argument 
based on the energy and transverse momentum conserva- 
tion. Assuming fii = — /i2 = /i > 0, a quasiparticle of en- 
ergy e, in either Gi or G2, may only have k = {iJ,+e)/fivF 
or k' = \fi — e\ /hvp as wavevector modulus. Conserva- 
tion of the transverse wavevector ky implies the Snell- 
Descartes law ky = k sin a = k' sin a' between the inci- 
dence angle a of the electrons and the reflection angle 
a' of the holes in Gi. Moreover a' is also the refrac- 
tion angle for transmitted electrons in G2. Since k' < k, 
choosing incident electrons with a above the critical angle 
ac(e) = arcsin(|(/i — e)\ /{fJ. + e)) yields a complete sup- 
pression of the Andreev reflection and electron transmis- 
sion [2^. Thus processes that are harmful for the CAR 
observation are both eliminated at once in channels with 
a > ac{s). In particular at e = /i, this suppression holds 
in all channels since adfJ.) = 0. Hence the whole current 
in G2 is purely carried by transmitted holes while the 
current in Gi is the superposition of the incoming and 



backscattered electronic currents. 

In order to investigate quantitatively the consequences 
of the previous Snell-Descartes argument, we consider a 
scattering state with an incoming electron in the conduc- 
tion band of Gi > 0) having energy e and transverse 
momentum ky. Owing to translational invariance along 
the interfaces, all scattered quasiparticle wavefunctions 
are expressed as 5'(a;)e*'^"^. 

We first consider channels with a below the critical 
angle ads) = arcsin(|(/i — £)| /{p + e)), or equivalently 
ky < k' . In the n-type graphene lead Gi, x < 0, the 
wavefunction is given by the following superposition of 
the incident electron, the reflected electron and the re- 
flected hole 



^-(2;) = (l,e^",0,0)e*'=™""^ 

-Kree(l,-e-*", 0,0)6-*'=™^"^ 



(1) 
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where Vee and The are respectively the amplitude for or- 
dinary and Andreev reflection at the Gi-S interface. The 
index a =sign(/^ — e) indicates whether the hole belongs 
to the conduction (cr = +) or the valence band (cr = — ). 

In the p-type lead G2 , a; > c?, the wavefunction consists 
in the superposition of the transmitted electron and hole 

■^{x) = tee(l,e-*''"',0,0)e*'''=''=°'"'(^-'') (2) 

where tee and the are respectively the amplitudes for 
elastic cotunneling and Andreev transmission (CAR) 
through the superconducting barrier. 

At incidence angles a > ads), namely for 
ky > k' , the expressions for the wavefunctions 
are still given by Eqs. (|l|2p except for the hole 
in Gi which is described by the evanescent wave 
rhe{0, 0, 1, iaOe'^'y'^'^'^ and for the electron in G2 de- 
scribed by tee{l, -icrC, 0, 0)e-V'=?-'^-"(^~''), where C = 
exp (arg cosh {ky/k'j). 

The wavefunction in the central superconducting bar- 
rier, < a; < rf, is the superposition of four kinds of 
waves given by a±,p(eT*'3, pe^F*'', e"*"*, pe-*'^)e''(*'=s±K)a;^ 

with p = ±1, ks — pLs/fivF 3> k,k' and k — 
■\/Ag — e^/Hvp- The phase (3 = arccos(e/Ao) is intrin- 
sically related to electron-hole conversion at a normal 
conductor-superconductor interface 

Demanding the continuity of the wavefunctions at x — 
and X = d yields the scattering amplitudes reef he, the, 
and tee (and a±^p) as functions of e,a,d and /i. In the 
limit d — > we recover the expressions for the transmis- 
sion and reflection amplitudes, tee and Tee, obtained so 
far in the study of the normal (non superconducting) n-p 
junction [3, [3], while rhe = the = 0. In the opposite 
limit c? 3> ^Oi the expressions for Andreev and normal 
reflection amplitudes rhe and Vee tend to those obtained 
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FIG. 2: Diagonal (top curve) and cross (bottom curve) dif- 
ferential conductances of the bipolar G1-S-G2 transistor as a 
function of the voltage for ^ — 0.3Ao and T = 0. 



in [9|, while transmission amplitudes are exponentially 
suppressed: t^e = the = 0. 

Diagonal and cross differential conductances of the 
G1SG2 heterojunction are deduced from an extended ver- 
sion of the Blonder-Tinkham-Klapwijk theory [2^. In 
the following, we assume that a positive bias Vi = V 
is applied to the normal lead Gi while the lead G2 and 
the superconductor S are grounded. Keeping in mind 
the critical angle effects discussed so far, the current li 
in the graphene lead Gi (i = 1,2) is represented as the 
sum of the currents if^ and !(' carried by channels with 
a < ac{V/e) and a > adV/e) respectively. 

We first obtain that the diagonal conductance dli / dV 
is finite at eV = fi for thin superconducting barriers 
d ^ ^Oi shown in Fig. [21 In contrast, the main charac- 
teristic of the GS contacts with infinite superconductor 
is the vanishing of the differential conductance at eV ~ [i 
i. 

We now consider the current I2 carried by electrons 
and holes transmitted in G2 when a positive bias is ap- 
plied to Gi. Channels with a < ads) contribute to the 
cross differential conductance as 



dV 



da 



k' 2 ' 2 

— \tee{s)\ COSa - \the{s)\ COS a ) , (3) 



where f ^ f{e ~ eVi) = l/(e(^-'=^i)/^ + 1) is the 
Fermi distribution of incident electrons in the lead Gi 
at temperature T. The factor 4 in = {Ae'^/h)N^ ac- 
counts for spin and valley-isospin degeneracy and Ng^ — 
{jjL -f e)W/{'KhvF) for a graphene sheet of width W . In 
contrast, the contribution to the cross conductance aris- 
ing from quasiparticles having a > ads) is always nega- 
tive 



''^^-fdsf-'^ 



dV 



ds 



Tr/2 



da\theis)\ cos a, (4) 



since then the electrons are evanescent waves which do 
not carry current. 




FIG. 3: Cross differential conductance of the bipolar G1SG2 
transistor as a function of the superconducting barrier width 
d at the "Fermi voltage" V — n/e and T — 0. The minima 
are located at d/fo = {n + l/2)nAo/fj,s- 



As shown in Fig. [51 the cross conductance 812/ dV ex- 
hibits a cusp at fj,/e being negative between Vd and Vc2 
and positive otherwise. This result may be understood 
further by comparing the cross differential conductances 
at Fermi bias eV = ^, at zero bias and at large bias 
eV » fj.. First at eV = m, dh/dV = dI>/dV is nega- 
tive for any width d because the critical angle vanishes. 
For voltages slightly shifted from fj,, the contribution 
remains dominant over owing to the larger angular 
integration interval in Eq.(|4]) compared to Eq.Q. On 
the contrary at zero bias, the critical angle is maximal, 
Q;e(0) = 7r/2, yielding dh/dV = dIf/dV. From the 
expressions of t;ie(0, a, d) and iee(0, a, d), one may show 
that the zero bias dl2 /dV is always positive. In conclu- 
sion, the cross differental conductance has at least a zero 
at a finite voltage Vd below eV = /i. A similar reversal of 
the cross conductance occurs at a voltage Vc2 above /i/e. 
The voltages Vd and Vc2 depend on the barrier width 
d , on /i and on ^s- 

The cross conductance is finite and oscillates as a func- 
tion of the superconductor size d as shown in Fig. [3l Re- 
markably, the lengths for which the conductance maxima 
occur are almost independent of fi. The experimental ob- 



servation of these oscillations requires Ac? <^ kg < d ^ 
^0 where Ac? is the typical fluctuation on d due to in- 
terface roughness. Owing to the good coupling between 
the superconductor and the atomic thick carbon layer, 
the Fermi wavelength kg^ is likely to be quite small in 
comparison to c? ^ ^0 01 ■ 

In addition, a recent experiment demonstrated that 
disorder ma y in duce spatial fluctuations of the chemical 
potential fj, 27|, |28[ . Since energy is still conserved, the 
general phenomena of AR and EC suppression at eV = fi 
should pertain although the wavefunctions are no longer 
plane waves. It should be very interesting to investigate 
the interplay of the AR and EC suppression with the 
formation of electron and hole puddles close to neutrality 
point [28| . 

Besides the intermediate energy regime fj, < Aq stud- 
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FIG. 4: Zero-temperature differential conductance dli/dV of 
the bipolar G1SG2 transistor in the hmits fi ^ Ao (left panel) 
and jj, — (right panel). The dashed lines are identical to the 
curves obtained in [|| for an infinite superconductor. 



led above, we now consider the extreme limits fi ^ 
and 11 = 0. Then the conductance dli/dV of a thin su- 
perconducting barrier (d ~ ^0) oscillates as a function of 
the bias voltage (Fig. 2]) due to the quasiparticles inter- 
ferences inside the superconducting barrier. In contrast 



are 



conductance oscillations in G1G2S junctions 
related to an interfacial barrier potential G2 separating 
Gi and S. Finally the cross conductance dlijdV is al- 
ways positive (EC dominated) because the phenomenon 
of EC suppression is lost at charge neutrality or when the 
Dirac points are largely outside the gap energy window. 

In conclusion, we have demonstrated that the favor- 
able kinematical conditions for splitting a Cooper pair 
towards two separate leads are met in a bipolar graphene 
transistor even in presence of weak disorder. This is the 
first step towards the realization of entangled states of 
massless electrons. Nevertheless clear-cut manifestation 
of entanglement depends on the actual relaxation and de- 
phasing mechanisms originating from intrinsitic effects in 
graphene as well as from the back action of the read-out 
devices. Finally, the proposed bipolar graphene transis- 
tor may serve as a very efficient Andreev beam splitter 
in Hanbury Brown-Twiss and Mach-Zender like experi- 
ments [3l| . 
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